Fluctuations in superconducting rings with two order parameters 
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Starting from the Ginzburg-Landau energy functional, we discuss how the presence of two order 
parameters and the couphng between them influence a superconducting ring in the fluctuative 
regime. Our method is exact, but requires numerical implementation. We also study approximations 
for which some analytic expressions can be obtained, and check their ranges of validity. We provide 
estimates for the temperature ranges where fluctuations are important, calculate the persistent 
current in MgB2 rings as a function of temperature and enclosed flux, and point out its additional 
dependence on the cross-section area of the ring. We find temperature regions in which fluctuations 
enhance the persistent currents and regions where they inhibit the persistent current. The presence 
of two order parameters that can fluctuate independently always leads to larger averages of the 
order parameters at Tc, but only for appropriate parameters this yields larger persistent current. In 
cases of very different material parameters for the two coupled condensates, the persistent current 
is inhibited. 

PACS numbers: 74.78.Na, 74.40.-n 



I. INTRODUCTION 

Fluctuations are extremely important near phase tran- 
sitions, and have therefore been the subject of intense 
research in the past. Particularly in superconductiv- 
ity, it has been shown that thermally driven electronic 
fluctuations, i.e. formation and dissociation of Cooper 
pairs close to the critical temperature Tc, can affect 
all relevant properties of a superconductor li Techniques 
that incorporate thermal fluctuations to the Ginzburg- 
Landau model are described in a recent review3 Fluc- 
tuations in mesoscopic loops are particularly interesting 
because their critical temperature is reduced in an oscil- 
latory fashion as a function of the magnetic field — a phe- 
nomenon known as the Little-Parks (LP) effect^ More 
importantly, as LP oscillations are directly related to flux 
(vorticity) entry in superconductors, one can identify the 
magnetic flelds for which fluctuations are particularly im- 
portant, as is the case of half-integer flux values.— The lat- 
ter experiment- detected current in the ring above Tc, a 
clear signature of fluctuations, and may be regarded as a 
paradigm for thin superconducting ring behavior, a case 
for which the theory for thermal fluctuations is known 
exactlyi^ The additional influence of quantum fluctua- 
tions was addressed in Ref . H. 

Superconductivity is essentially a macroscopic quan- 
tum state with long-range phase coherence, therefore de- 
scribed as a single wave function. Superconductors with 
several order parameters have recently attracted great 
attention due to the discovery of MgB2 and high Tc su- 
perconductivity in pnictides. In such cases, thermal exci- 
tation allows contributions from multiple wave functions 
and one may expect a dramatically different behavior of 
the system. With that as motivation, we here explore the 
interplay of the wave functions and thermal fluctuations 
in superconducting rings with two order parameters. The 
rings we will consider need not be made of a two-band 



superconductor, but may also consist of two thin super- 
imposed superconducting rings possibly separated by 
an isolating layer, such as the active part in readily made 
experiments with annular Josephson junctionsi^ 

We conduct our theoretical analysis in the framework 
of the Ginzburg-Landau (GL) theory. The multiband 
GL equations were developed long agoj^ in the case of 
two bands the free energy density has the form 



i/=l,2 \ 



2 + ||A.|4 + 7^,|nA.|2 



-7(AiA* + A2AI), (1) 



where Ai^2 are the order parameters, ai.2, bi^2, ^1,2 
and 7 are material parameters and 11 = V -I- 27riA/$0j 
with A the vector potential and $0 the superconduct- 
ing flux-quantum. Zhitomirsky and Dac^*^ obtained ex- 
pressions for the material parameters in a multiband 
superconductor using Gor'kov's technique. Kogan and 
Schmalian^^ recently emphasized that consistency im- 
poses conditions on the temperature dependence of these 
coefficients, which results in the same coherence length 
for both order parameters in a two-band superconduc- 
tor. Shanenko et ali^ went on to show the importance 
of terms of higher order in temperature, and the resulting 
separation of characteristic lengths for the two bands. We 
should note, however, that fluctuations move the order 
parameters astray from equilibrium, so that in general 
their ratio is not constant. Moreover, besides two-band 
superconductors, we are interested in relating our results 
to additional systems, where the coefficients in Eq. ([T]) 
may have a different temperature dependence. In this 
paper we thus adopt the standard GL approach, where 
the material parameters in Eq. ([I} are arbitrary functions 
of the temperature and any required restriction will be a 
particular case. 
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II. METHOD 

In this paper we deal with one dimensional supercon- 
ducting rings with two order parameters, extending the 
results obtained by von Oppen and Riedel^ from single 
to two order parameters. Instead of the formalism of 
Ref.H, we will follow a slightly different approach, which 
in our view is conceptually simpler. We start by absorb- 
ing the coefficients Ki^2 into the order parameters and by 
switching to a gauge invariant formulation, i.e. we define 



A^{0) = exp 



/ 2T:iR 



A{e')de' \^K,A,{e) , (2) 



where R is the radius of the ring, 9 is the angle along the 
ring and A is the tangential component of A. Likewise 
we define 



0,1, , K 



1 



(3) 



With these definitions the free energy density becomes 



/ 



1^=1.2 



aJAJ 



1 



d9 



-7(AiA; -t- A^A^) . (4) 



— 1/2 

With the normalizations we are using, ai, is the coher- 
ence length for the order parameter A^ in the absence of 
coupling to the other order parameter. As in Ref. @, we 
consider a uniform ID ring, which in particular does not 
alter the applied field (screening is negligible), so that no 
magnetic energy has to be added to /. 

In order to have a more intuitive picture of the prob- 
lem, we represent the complex order parameters Ai, by 
two-dimensional real vectors r^, such that in polar coor- 
dinates r^, = \A^\ and the angular coordinate -d^, is the 
phase of A^. Integrating / over the volume of the ring, 
its free energy becomes 
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where w is the cross section of the ring and 

R^ f 9 9 . ^1 4 b' 



(5) 



V 



2 ' 2 
2^rir2 cos(i?i - ih) 



(6) 



As in Ref. 5, ri,2(^) may be regarded as the trajectories 
of two fictitious particles during a period of time — tt < 
6 < TT. The first two terms in the integrand of Eq. ([5]) 
then represent their kinetic energy and V their potential 
energy. 

Following Ref. [l^, a pair of functions r^{9) is inter- 
preted as a microstate of the system and F as the energy 



of the system for that microstate. It follows that up to an 
irrelevant multiplicative constant the partition function 
is 



J VyiVv2 exp{~F/kBT) , 



(7) 



where T is the temperature and / X'ri2?r2 denotes inte- 
gration over all functions (9) with appropriate period- 
icity. Since A^ are single valued, r,^{9 = tt) = r,^{9 = — tt) 
and = tt) ^ d^{9 — — tt) + 2ttlp, where ip is the flux 
enclosed by the ring divided by ho- 
using slight adaptations of Eqs. (2.14), (2.15), (2.16), 
(2.22) and (2.23) in Ref. [H, Z can be brought to the 
form 

Z= J2 exp(-27re„/^) 

n 

X J dridr2*:[r,(0 ^)]*„K(0 = -tt)] . (8) 

Here e„ and [r^] are a complete set of eigenvalues and 
normalized eigenfunctions of the fictitious Hamiltonian 



(9) 



where the Laplacian acts on r^, and S = kBTR/2w. S 
has dimensions of surface in the plane of the trajectories 
of ri^2 and dimensions of force in reality. The integral 
in Eq. is taken over the entire planes of motion for 
each particle, but for every argument Vi, in vlf„ we have 
to take the corresponding argument in ^* . 

We note now that the angular momentum opera- 
tor Lz — —i{d/ddi -I- 8/8^2) commutes with H. We 
can therefore choose the set of eigenfunctions {v^n} 
with well defined angular momentum, i.e. they can 
obey Lz'ifnj — i'^n.i and therefore have the form 
^'n,£(ri,r2) = Y^ti '^n,i,iA'>'i,r2,'&i-'d2)exp[i{eidi + {e- 
(^1)1^2)], with £ and £1 integers. In view of the peri- 
odicity conditions of r^, it follows 'i>n^i[r^{9 = tt)] — 
exp{2Tri£ip)'i'n.e[i^i^{9 ~ — tt)]. By substitution of this re- 
sult into Eq. ([SJ we obtain 



with 



Z = cxp{—2Tri£ip)Zi 



Zg = y^exp(-27r£„^f/6') , 



(10) 



(11) 



where summation in Eq. (llOp is made over all integers 
and in Eq. ITTI) over all the states with total angular 
momentum £. Since H is invariant under the transfor- 
mation {'i?i,i32} {— "i?!, — ^^2}, £n.-i = Sn,i and we can 
also write 



Z = Zo + 2 ^ cos{2Ti£ip)Zi 



(12) 
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Once the partition function is known, all the equilib- 
rium quantities can be derived from it. The average cur- 
rent around the ring is 



{I) 



kBTd\nZ_ keTin 
$0 dip $0 Z 



sm{2TTlip)£Zi . (13) 



Following similar steps to those that led to Eq. (IT^ . we 
obtain that the statistical average of any function of the 
absolute values of the order parameters is 



(p(lAi|,|A2|)) = i 



1=1 



(14) 



where Pi = X]„ exp(-27re„,£/S')(nf |p(ri, r2)In£). Here 
we introduced the matrix element {nt\p{ri,r2)\'nl) — 
J dridr2p(ri, r2)|4'ra,^(ry)p, which may in practice be 
evaluated in any convenient basis. 



(iv) the total angular momentum equals the sum of the 
angular momenta of each particle and the total energy 
equals the sum of their energies. As a consequence of 
these features, the partition function in Eq. (ITUl) becomes 
Z = where Z^"'^ is the value of Z obtained when 

is replaced with the value that corresponds to parti- 
cle v only. It follows that equilibrium quantities such as 
the average current or the average energy will equal the 
sum of the separate contributions of particles 3 and 4. 

The Hamiltonian of the fictitious particles 3 and 4 is 
just that of two decoupled harmonic oscillators and its 
eigenvalues are well known: elf f = (S'i?/e3.4)(2n + |£| + 
1). The sum in Eq. (jlip becomes a geometric series and 
we obtain 



.(3,4) _ 



,-27r|£|_R/C3,4 



2sinh27ri?/^: 



3,4 



(18) 



Substitution into Eq. ([TU]) gives 



III. EVALUATIONS 



A. High temperature 



Let T be sufhciently higher than Tc, so that the order 
parameters are small and the quartic terms in V can be 
neglected. For high T we can also assume ai + a2 > 
a/ (oi — 02)^ -I- 47^ and define the quantities 



V = 



ai - 02 



?3,4 



2 - «2)2 + 472 ' 
2 



ai + a2 T - 02)^ +472 



(15) 



— 1/2 

For 7 = and ai > a2, ^3.4 = 12 1 i therefore re- 
gard ^3_4 as a sort of coherence lengths in the presence of 
coupling. 

We can now define a rotation in the 4D space of both 
particles through 



1 



r2 = 1-3 



1 

2 + ^ I"4 , 



(16) 



With this transformation, 
Eq. ^ becomes 



the "potential energy" in 



quad 



^{R'/2)[{rs/^sr + {n/^,Y 



(17) 



The following features should be noted: (i) in the co- 
ordinates {r3,r4}, Vf -I- V2 still has the meaning of the 
Laplacian in the 4D space; (ii) since on passing from 
6 = —IT to 6 — TT both and '§2 increase by 2'Kip, and 
since r^^^ are linear combinations of ri_2 with fixed coef- 
ficients, also their angles ^3^4 increase by 2Trip; (iii) the 
angular momentum operators —id/d'd^ ^ for each sepa- 
rate particle now both commute with the Hamiltonian; 



Z(3'4) ^ [2(cosh27ri?/,^3,4 - cos27rv3)]"i (19) 
and the average current in the ring equals 



(-^quad) — 

-(27rsin27r(/5fcBT/$o) [(cosh27ri?/^3 
+ (cosh 27ri?/^4 — cos 27r(^)^"'"] . 



cos 2'KLp) ^ 
(20) 



From Eq. pO|) we can obtain the Little-Parks temper- 
ature, i.e., the temperature for the onset of supercon- 
ductivity in the absence of fluctuations. Without fluc- 
tuations the current vanishes above this onset; this is 
implemented by taking the limit ksT — > in the first 
factor in Eq. (|20p . At the LP temperature the current 
becomes nonzero, requiring divergence of the second fac- 
tor, i.e. ii?/^3 — if mod 1. From here, the LP condition 
is 



R'Wiai - a2f + 472 - (ai + 02)] = 2^"" , (21) 

where for simplicity of notation we restrict ourselves to 
the range |(^| < 1/2. Near we can write ai, = a^c—OivT, 
with r = (Tc — T)lTc and aica2c = 7^; if in addition 
{if/RY < "1,2, condition ([H]) is fulfilled for 



TLP 



flic 



a2c 



aicCt2 + a2cQ!i i?2 



(22) 



B. Hartree approximation 



In the potential given by Eq. ([6]) we now make the re- 
placement (6i.2/2)rf 2 bi_2{r\ 2)'^! 2- high tempera- 
tures both terms are negligible and at low temperatures, 
where fluctuations can be neglected, they both lead to 
the same "force" -WM- 

With this approximation the potential becomes again 
quadratic, so that we can still use the results of the pre- 
vious section by substituting 77 — 77' and ^3^4 — >■ ^34, 



4 



with 



M _ 

S3, 4 — 



2yK" 

2 



472 



ai.2 = ai,2 +^1.2(^1.2) . 



(23) 



In order to implement this approximation, we have to 
evaluate {r\^. {r\ are given by 



'3,4/ 



^^3,4 ainZ(3-4) 

~2^9(T7iiJ 

sinh27ri?/^^ 
R cosh 27r_R/^3 4 — cos 27r(/? 



(24) 



on the other hand, since (ra • r4) = 0, {rX^j = (1/2 =F 
r?')(r|) + (l/2±77')(r|), hence {rl^) = (l/2Tl/4r,')(r?) + 
(1/2 ± l/4r/')(?'i)- Substituting this into Eq. ^ leads 
to 



^^47) ^'^^^+G±47 



3,4 



sinh27ri?/^^4 



i? cosh ^-KRji^^ 4 — cos 27r(/j 



(25) 



This is a system of two equations for obtaining (rf ) and 
krl). 

For a general situation, Eqs. ([25]) have to be solved 
numerically, but we can find asymptotic expressions for 
some special situations of interest. Sufficiently above the 
critical temperature one usually has 2itR/^'^ 4 ;» 1, so 
that the fractions at the right of Eqs. reduce to 1. 
Near Tc and for a range in which hv{r1) ^ ai/|r| ^ 
we obtain 



(|Ai,2p(T)) = (r?,2) 
keT ( 



2(aic + a2c)^/^w 



flic, 2c ■ 



a2c,ic(aic + a2c 
\J {a\ca.2 + a2cai)|T 




It is interesting to note that (|Ai^2p(7')) decreases very 
moderately with T — Tc. If the radius of the ring is suf- 
ficiently large, we still have 27r_R/^3 4 1 at T = T^, 
assuming that b„{r1) ^ stiU holds leads to 



(|Ai,2r(rc)) ^ 

aic,2c[3/3i,2 



^2 T 



1/3 



I 4(aic + a2c)/3i,2w2 
a2c,ic(^i + b2)]kBTc 



6(aic + a2c)^/^/3i,2W 



(27) 



with /3i,2 = &i,2a2c,ic + fe2,iaic,2c/7^- Comparison of 
Eqs. (gni) and (gT]) with the values of |Ai_2|^ in the ab- 
sence of fiuctuations enables us to estimate the range of 
temperatures for which fluctuations are important. 
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FIG. 1: Average Cooper-pair densities | A„=i.2(=a-,7r) |^ in the 
ring as a function of temperature, calculated in the Hartree 
approximation. Each |Ai/p has been multiplied by Pi, for 
the purpose of comparison with the values in the absence of 
thermal fluctuations (dashed line). The sample is a ring of 
radius 10^* cm and cross section 10~^° cm^ that encloses no 
magnetic flux. The material parameters are those of MgB2, 
taken from Ref.[il. Inset: {R/^'sf for the same rmg. 



Figure [T] shows the values of (A^ 2) = (r^ 2) near Tc for 
MgB2, using the material parameters reported in Ref . [T^. 
In the absence of thermal fluctuations, /3i,2(A4 2) is given 
by the dashed line for both bands. The inset in the figure 



shows the behavior of R/^'^ near Tc 
practically constant in this range. 



?4 ?3 remains 



Exact evaluation 



Y(V? + V^) + ^(fc?r? + fc2S'),(28) 



We define the basis Hamiltonian 
HB{ki,k2) - 
which has the known set of eigenfunctions 

^„„.„„.,..(ri,r2) = C n rl^'^le-^'^-^/^^ 

u=l,2 

X iF^{-n,,\l,\ + \,RkyjS)e''-^^ (29) 

with eigenvalues i?S' [fci (2ni -h |^i | 1) (2n2 |-^2 1 + 1)] • 
Here C is the normalization constant and iFi is Kum- 
mer's hypergeometric function. We can then evaluate 
any matrix element Hi j = {tpi\H\7pj) , where the func- 
tions '0i=i,...,7v are the functions with lowest eigenvalues 
within the basis of the Hilbert space provided by Eq. (|29p . 
More precisely, when evaluating Z^, we include in the set 
{iJi] only eigenfunctions that obey £1+^2 = ^- If N is 
sufficiently large, then the lowest eigenvalues of the ma- 
trix (Hij) will be a sufficiently accurate approximation 
for the lowest eigenvalues of the operator H. In practice, 
rather than fixing N, we fix maximum values for rti, \£i\, 
n2 and |^2|- 

We are interested in choosing ki and k2 so that an 
accurate approximation is obtained without N becoming 



-0.010 -0.005 0.000 0.005 0.010 0.015 0.020 
(T-T,)IT, 

FIG. 2: Average current in a two-order-parameter supercon- 
ducting ring as a function of temperature. For the top curve 
the parameters are as in Fig. [T] except that the normalized 
flux is </5 = 0.25. For the lower curves the cross section w of 
the ring is smaller. The current was evaluated using a trun- 
cated basis of the Hilbert space, provided by the functions 
in Eq. (|29|) with maximum quantum numbers ni.max = 11, 

1-^1 jmax = 22, n2,max = 4, and |^2lmax = 8. 



prohibitively large. From Eq. (|29|) we see that the scale of 
{r1) is given by S/Rk^; we therefore set ki, = pS/R{rl), 
where (r^) is obtained from the Hartree approximation 
and p is still a free parameter. Since for a good approx- 
imation the eigenvalues should actually be independent 
of p, we mimic this situation by minimizing the lowest 
eigenvalue of (Hij) with respect to p. 

Figure [5] shows the currents as functions of the tem- 
perature obtained with this method for rings of three 
different cross sections, with the material parameters of 
MgB2. For completeness, the figure includes also values 
of current that are too small to be experimentally ob- 
servable. The temperature range covered here is much 
wider than the range presented in Ref. [^, where the tem- 
perature scale is given by the Thouless correlation energy 
(divided by ks) which is of the order of the LP tempera- 
ture; for a MgB2 ring with radius of the order of a micron 
the LP temperature is of the order of 10~^ Tc- 

For sufficiently high temperature (the required tem- 
perature decreases with cross section w), the current be- 
comes independent of the cross section. For the parame- 
ters taken here and T = 3Tc, a ring with w = 10~"'^''cm^ 
and a ring with w — lO^^cm^ carry the same average 
current (within 1% difference). At the other extreme, 
far below Tc, i.e. where fluctuations are unimportant, 
(/) oc w. However, the dependence of (/} on w near 
Tc is not intermediate: we notice in Fig. [2] that decrease 
of w by an order of magnitude near Tc leads to a current 
decrease of several orders of magnitude, whereas for in- 
termediate behavior the current would decrease by one 
order of magnitude at most. Figure [3] shows the scaled 
current against the scaled temperature for the same rings 
as in Fig. [51 close to Tc- We notice that, in spite of the 
moderate influence of temperature on (|Ai.2p} predicted 



-1.0 -0.5 0.0 0.5 1.0 1.5 

10V-''(7'-7'..)/7-, (cm-/3) 

FIG. 3: Scaling of the function {I{T)) with the cross section. 
The parameters are the same as in Fig.[2l lo = 10"^" cm^; 
Q: w ^ 10"" cm^; D: w = 10"^^ cm^ The straight line is a 
guide for the eye. 



by Eq. ([5^ . the current decreases exponentially. We also 
find that for smaller cross sections the rate of change of 
(/) is slower. We empirically found that the scaling w^^"^ 
leads to a universal curve. 

We attribute didactic interest to understanding the be- 
havior of our methods far below Tc- There, convergence 
of the series in Eqs. (fT2l) and (fT3)) becomes slow, and 
numeric implementation of the exact evaluation becomes 
inefficient. Below certain temperature, in Eq. ([55]) 
becomes negative, and interpretation of the potential in 
Eq. (flT)) as that of a harmonic oscillator, and the sum 
of convergent geometric series that led to Eqs. ([50)) and 
([55)1 is no longer justified. Nevertheless, the expressions 
in Eqs. ([50)) and ([55)) are analytic functions of ^3^, than 
remain meaningful and are expected to remain valid be- 
yond the range in which they were proven. Numeric im- 
plementation of the Hartree approximation requires spe- 
cial care in order to pick the relevant rather than spurious 
solutions of Eqs. ([55)) . In the limit T ^ 0, R/S,'^ i(p, so 
that the right hand side in the first of Eqs. ([55)) does not 
vanish. Taking the appropriate limits in Eqs. ([50)) and 
()25p we obtain that the current in the Hartree approxi- 
mation is Ih{0) — —4:7Tw{Al + A2)f/R^o, exactly as in 
the absence of fluctuations. 

We conclude this section with a review of the accu- 
racy of our evaluations. The accuracy of the "exact" 
evaluation can be estimated by repeating it with reduced 
maximum values for and \£^\. We found the largest 
inaccuracy for large w and small T. In the results pre- 
sented in Fig. [51 the maximal inaccuracies are of the order 
of 10%. Figure [4] compares our approximation methods 
against the exact evaluation for w = 10~^'^cm^. One can 
see in the flgure that all the approximation methods are 
very inaccurate precisely in the most interesting region, 
i.e., close to Tc- The range of temperatures where the 
approximations are inaccurate is larger for smaller cross- 
sections w. Note also that there exists a range of tem- 
peratures for which the mean fleld current is larger than 
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(r-r,)/r,. 

FIG. 4: Ratio of the current values obtained with approxi- 
mated methods to the exact evaluation. Parameters are the 
same as in Fig. [TJ except that here ip = 0.25. G: Gaussian 
(quadratic) approximation; H: Hartree approximation; MF: 
mean field evaluation (i.e., without fluctuations). The de- 
scent of the MF curve to at the Little-Parks temperature 
looks vertical in this scale. 
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FIG. 5: Comparison among currents as functions of the tem- 
perature for systems with different coupling 7, but with oth- 
erwise identical parameters. The value of 7 is marked next to 
every curve in E-notation (e.g. 3E7 denotes 7 = 3 x lO^cm"^). 
The other parameters are R = 10~^ cm, w = 10"^'' cm^, 
a = lO^^cm-2, b = 10^'' erg" ^ cm" \ ksTc = 10"^^ erg, 
= 0.25. 



our exact evaluation, meaning that thermal fluctuations 
inhibit the current in this region. 



IV. INDEPENDENCE AND ASYMMETRY OF 
THE ORDER PARAMETERS 

The most conspicuous qualitative differences of a two 
order parameter system, as compared to a system with 
a single order parameter, are independence and asym- 
metry. By independence we mean that at a given point 
and time the two order parameters are not necessarily 
equal to each other; by asymmetry we mean that the av- 
erage values of the order parameters are not necessarily 
the same. In this section we investigate the influence of 
these properties. 



A. Symmetric case 

We start by considering independence while assuming 
equal coefficients for both order parameters. For simplic- 
ity, we assume 

fly = 7 - ar , , (30) 

with 7, a and b constants. As an illustration, we 
may think of a film of a uniform single-parameter ma- 
terial of thickness zq with energy density — 2aT|Ap -|- 
6|A|* -I- 2|VA|^. If we decide to denote by Ai (resp. 
A2) the value of A in the upper (resp. lower) half of 
the film, substitute the z-derivative by a finite differ- 
ence and average over zo, we obtain the energy density 
-ardAip + lAaH + b{\^,\^ + |A2n/2 + |V,,Ai|2 + 
|V,,A2|2 + (8/z2)(|Ai|2 + IA2I2 - AiA* - A2At), with 



^ xy being the component of the gradient in the plane of 
the film. One can easily identify that we have recovered 
the energy density for two order parameters, with cou- 
pling 7 = 8/zq. In the limit zq ^> 0, 7 — > 00, and Ai and 
A2 are the same; in the opposite extreme, 7 — >■ 0, and Ai 
and A2 are independent, while in the general case they 
are correlated. Following this analogy, it is very easy to 
obtain results for the cases 7-^-0 and 7 00: the case 
7 is equivalent to that of two single parameter sys- 
tems in parallel, and the case 7 00 is equivalent to 
that of a single parameter system with a doubled cross 
section. 

Figure [5] compares calculated average currents {I{T)) 
as the parameter 7 is varied in the range < 7 < 00, 
while all the other parameters are common to all curves. 
For a facilitated comparison, all the functions have been 
divided by {Ioo{T)), the current obtained for 7 — > 00. 
Figure [5] shows the temperature range close to T^, far 
below Tc the influence of fluctuations is negligible and 
all the curves should coalesce. For every temperature, 
we note that as 7 increases from to 00, {I{T)) changes 
from (/o(T)) to {Iao{T)). However, this change is not 
monotonic: |(/(r))| initially decreases and after reaching 
a minimum increases towards |(/oo(7^))|. The fact that 
|(/o(r))| < |(/oo(T))| for T w Te may look surprising, 
since 7 = means larger freedom than 7 — > cxi and we 
would therefore expect larger fluctuations in the former 
case. We will see in the following that indeed the order 
parameters assume larger values for 7 = 0; however, they 
may be less coordinated, resulting in a smaller current. 

The solid curves in Fig. [5] show the values of (p for which 
the current is maximum aX T — Tc for the limiting cases 
7 = and 7 — > 00. From a dimensional analysis we find 
that in the present situation the temperature enters the 
operator H/ S only through the combination kBTjDR^ /w, 
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FIG. 6: Magnetic flux for which the fluctuation current is 
maximal, and value of the order parameters at T = Tc, as 
functions of T^. (|Ap) wcls evaluated at — v^max- Each 
curve is marked by its value of 7 and by an arrow that points 
to the relevant j/-axis. The thin straight line (red online) 
highlights the asymptotic power dependence of bR^{\Af) on 
kBTcbR^/w. 



SO that ipmax is a function of this quantity. Since the 
mean field current has its maximum at (pmax > 1/4, it 
is interesting to note that ipynax can be smaller than 1/4. 
The curve for 7 — > 00 can be inferred from the case 7 = 0: 
in order to obtain it at a given Tc, we have to double 
the value of w. Since (pmax is a function of kBTcbR^/w, 
doubling w is the same as dividing Tc by 2, i.e., at a 
given Tc, the value of ipma.x for 7 — > cxd is the same that 
fmax for 7 = had at half that temperature. With a 
logarithmic x-axis, this relation gives a shift of the curve 
to the right. 

The dashed fines in Fig.lHlsfiow the values of 6i?^(|Ap) 
at T = Tc and — ipma.x, evaluated by means of 
Eq. (HH). Except for ksTc < w/bR^ or keTc > w/bR^, 
we obtain (|A|2) « Qm[bR^)-^{kBTcbR^/wf/^ = 
0.68(fc|T2/6w2)i/3 for 7 = 0. For 7 ^ 00 (|A|2) is 

smaller by a factor 2"^/^. The first term in the Hartree 
approximation value in Eq. (|27l) is smaller than the result 
obtained for 7 00 by about 7%. 

Figure [7] shows the average current evaluated at T = 
Tc and = v^max- As already found in Ref. {I{Tc)) 
is not a monotonic function of Tc, but has a maximum 
instead. As discussed above, the curve for 7 ^> cxd is 
obtained as a shift of the curve for 7 = 0. What we 
learn from this curve is that for ksTc < 0.163w/bR'^ 
\{Io{Tc))\ > |(/oo(T'c))|, meaning that independence of 
the order parameters enhances the fluctuation current, 
whereas the opposite occurs for fc^Tc > 0.163w/bR^. 



B. Asymmetric case 

There are three material parameters that can differ be- 
tween the order parameters: aic ^ a2c, O-i 7^ a2 and &i ^ 
62. Since near Tc the usual case is |ai,2'''|, &i,2|Ai^2p ^ 
|aic — o-2c\, we focus on the influence of the difference 




0.00 t , , , , , L 

0.01 0.03 0.1 0.3 1 3 



FIG. 7: Maximum current at T = Tc as a function of the 
scaled Tc for 7 = and for 7 — )■ cxd. 
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FIG. 8: Ratio between the fluctuation currents for the case 
tsic 7^ a2c (denoted I^) and the case aic = a2c (denoted 7=). 
For the solid lines b = 3x lO^^^erg-^cm-^ {ksTcbR^ /w = 0.03) 
and 7 = 3 X 10®cm~^; for the dashed lines b — 10^^erg~^cm~^ 
{ksTcbR^ /w = 1) and 7 = 3 x lO^cm"^. In all cases <p = 
</5max. The other parameters are the same as in Fig. [S] 

between aic and a2c- 

Figure [8] shows the ratio between the fluctuation cur- 
rents for the cases aic 7^ 0,2c and aic = o-2c, while all 
the other parameters are kept unchanged. Although the 
values of kBTcbR^ /w and 7 do have some influence, the 
general trend is that the difference between aic and a2c 
inhibits fluctuation supercurrent in the region T k, Tc, 
with this effect being stronger for T < Tc- 

V. CONCLUSIONS 

Motivated by recent surge in interest in the physics 
of coupled condensates in two-band superconductors, we 
have analyzed the role and importance of fluctuations 
in superconducting rings with two order parameters. We 
have extended the analysis of the fluctuative regime made 
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by von Oppen and Riedel^ for the single order parameter 
case, based on the Ginzburg-Landau energy functional. 
Further, we have made semi-analytic evaluations of the 
influence of fluctuations on the persistent current and on 
the order parameters in the ring, as functions of temper- 
ature, coupling between the order parameters, and mag- 
netic flux. We have identified the ranges of parameters 
where fluctuations inhibit or enhance the persistent cur- 
rent in the ring, and pointed out the influence of the cross 
section of the ring as well as the influence of the freedom 
of the order parameters to undergo separate fluctuations. 
Although the influence of fluctuations is most important 
close to Tc, we have also studied the behavior far from 
Tc, providing a complete picture. 

In addition to two-band materials, our findings apply 
to artificially made systems of two superimposed rings, as 
encountered in experiments that involve annular Joseph- 



son junctions. The present study can also serve as a 
guideline for theoretical efforts and interpretations of ex- 
perimental data in systems described by multiple order 
parameters in the fiuctuative regime, especially in nan- 
othin samples, which are always effectively multiband 
due to quantum confinement^^ and where fluctuations 
are of outmost importancej^^ 
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